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Abstract. We present a derivation of the Gribov equation for the gluon/photon Green's function D{q). Our 
derivation is based on the second derivative of the gauge-invariant quantity Trln_D(g), which we interpret 
as the gauge-boson 'self-loop'. By considering the higher-order corrections to this quantity, we are able to 
obtain a Gribov equation which sums the logarithmically enhanced corrections. By solving this equation, 
we obtain the non-perturbative running coupling in both QCD and QED. In the case of QCD, as has 
a singularity in the space-like region corresponding to super-criticality which is argued to be resolved in 
Gribov's light-quark confinement scenario. For the QED coupling in the UV limit, we obtain a oc 
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behaviour for space-like = 
theory in the UV limit. 



-q . This implies the decoupling of the photon and an NJLVL-type effective 



PACS. ll.lO.Lm Nonlinear or nonlocal theories and models - 11. 15. Ex Spontaneous breaking of gauge 
symmetries - 11.15.Tk Other nonperturbative techniques - 12. 38. Aw General properties of QCD 



1 Introduction 

Gribov's programme [ll[2] for dealing with the problem 
of strongly interacting quarks is based on the picture of 
super-critical binding of light quarks, leading to chiral- 
symmetry breaking and confinement through the reorga- 
nization of the quark (the Dirac-sea) states in the vacuum 

m- 

The central analytical tool used in this approach is the 
Gribov equation. 

In order to derive the Gribov equation, let us con- 
sider the (all-order) self-energy correction, S{q), to the 
quark propagator G{q). We denote the 4-momentum, g^, 
derivative d/dq^j, by 9^. The raising and lowering of in- 
dices is implicit. The 4-dimensional Laplacian operator is 
then defined as = 9^9;^, using the Einstein summation 
convention. 

The application of to the all-order S{q) yields an in- 
finite series of perturbative diagrams, in which some quark 
or gluon propagator is replaced by its derivative. 

The leading, logarithmically-enhanced, contribution to 
S{q) comes from phase-space regions with large hierar- 
chies of internal momenta. Corresponding to this region, 
the leading contribution to the second derivative of S{q) 
comes from diagrams in which both derivatives act on the 
same propagator, i.e., d^D{q — q') where D{q — q') is the 
internal gluon (in this case, bare) propagator carrying 4- 
momentum q—q'. Neglecting the vacuum polarization con- 
tributions, it is possible to arrange the internal momenta 
in such a way that derivatives of the quark propagator 
G do not occur. This simplifies things (in the Feynman 



gauge) since we have the relation: 



{q — q'Y + ie 



(1) 



which removes a loop integration and make it formally 
equivalent to the amplitude for the emission of two zero- 
momentum gluons. 

The leading contribution to this double emission is ob- 
viously /)jG((7)/)j, i.e., the successive emission of two glu- 
ons from a quark line. The 0-momentum emission vertex 
is related to the Green's function by the Ward-Takahashi 
identity: 

r^{q,q,Q)^-d^G'\q). (2) 

Since S{q) gives the running of G~^{q), the end result 
of this discussion is the following Gribov equation which 
sums all the leading logarithmically enhanced contribu- 
tions: 

d^G-\q) = gd^G-\q)G{q)d^G-\q) + 0{g^). (3) 

The coupling g is a/ it for QED and CFCts/T^ for QCD. 
Because of the removal of the loop integration, this equa- 
tion is local in the momentum space, and benefits from 
the absence of the momentum space integration over the 
dangerous IR (infra-red) region. 

As mentioned above, the renormalization of the gluon, 
or the photon in QED, propagator has been neglected in 
this discussion, though it is possible to incorporate this 
effect partly by replacing g in eqn.[3]by a running coupling, 
.9(9). 
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Leaving aside the problem of the coupled evolution 
of the quark and gluon Green's functions for now, the 
problem is to write down an analogous Gribov equation 
for the gluon/photon sector. 

In this case, the vertex r{q,q',q — q') is evidently 
no longer constant at the tree level, and is linear in the 
momenta. Therefore the simple application of the above 
method, of applying to the vacuum-polarization opera- 
tor Ilxa-, is insufhcient in the sense that it yields extra con- 
tributions due to the derivatives of the vertex. Gribov's so- 
lution to this problem [Ij was to employ the linear Duffin- 
Kemmer formalism, in which the interaction is constant, 
so that we can analyze the problem in a similar way to 
that of the fermion sector. However, the resulting equa- 
tion, being a coupled second-order differential equation, is 
hardly manageable, and furthermore suffers from artificial 
divergences which plague the Duffin-Kemmer formalism 

On the other hand, the analysis based on the Duffin- 
Kemmer formalism is illuminating at least in the sense 
that it illustrates how to recover the Ward-Takahashi iden- 
tity, eqn. [21 in the zero- momentum limit. This problem 
and the problem of gauge fixing are related. These com- 
prise the main difficulty in formulating a Gribov equation 
for the gluon/photon sector. 

On the other hand, for the problem of the UV (ultra- 
violet) evolution of the QED coupling, Gribov postulated 
[ll|2l[4] a formulation based on the third derivative of the 
polarization operator 11 Three differentiations are nec- 
essary to remove the UV divergence. In this case, Gribov 
obtained: 




^ = In I I , and Ff is the fermion anomalous dimension 
dlnZf /d^. In the weak-coupling limit, we have Fq = 
and so the right-hand side gives —2/3. The ordinary RGE 
(renormalization-group equation) evolution is recovered in 
this case. In the strong-coupling limit, Ff tends to (1 + 
•\/3)/2, and so the right-hand side tends to zero. In this 
case, the running coupling continues to evolve, but more 
slowly than in the RGE evolution which suffers from the 
problem of the UV Landau pole. 

In the UV limit, the right-hand side of eqn. 2] tends to 
—2/g'^. This leads to a logarithmic evolution, g ^ ^. Hence 
the QED (actually U{1)y) couphng grows large but finite 
for finite ^ and, as postulated by Gribov [5], causes the 
formation of super-critical states, which could be identified 
with the Higgs and Goldstone bosons in EWSB (electro- 
weak symmetry breaking). 

In the derivation of both eqn. [3] and 31 the running- 
coupling effects due to the renormalization of the internal 
boson lines are neglected. Furthermore, in the QED case, 
the right-hand side of eqn. [H is a derivative of a phase- 
space integral, and so the formulation is not strictly local. 

The purpose of this paper is to present an approach 
to deal with these problems. Starting from a gauge invari- 
ant expression which corresponds to the gauge-boson 'self- 
loop', and collecting the logarithmically-enhanced terms 



as in the derivation of eqn. [31 we obtain a Gribov equa- 
tion for gauge bosons, and an evolution equation for the 
running coupling, in a closed form. 

It is found that the equation can be integrated analyt- 
ically. We analyze the solution for both QCD and QED. 

In QCD, the solution exhibits a branch-cut singularity 
at space-like momentum = — 9^ = ^qcD' indicating 
vacuum instability due to the formation of the critical 
state. The problem of vacuum instability is resolved in 
Gribov's scenario due to the reorganization of the quark 
states which is best understood using the Dirac-sea pic- 
ture. We argue that this will remove the singularity and 
give rise to a Green's function for the gluon which has 
the analytical properties that are expected for confined 
particles. 

If, on the other hand, there are no quark states which 
become super-critical, there is a problem as to how to sta- 
bilize the vacuum, and we believe that this is not possible 
within the framework of ordinary field theory. 

In the case of QED at high energy, we found that the 
coupling grows linearly with = —q^. This corresponds 
to the decoupling of the photon, and an effective high- 
energy theory described by a contact interaction term, as 
in the NJLVL model ^6]. This evolution of the coupling 
is faster than the logarithmic behaviour found by Gribov. 
We show that provided that one includes the contribution 
due to internal photon renormalization in eqn. [4l the oc 
behaviour of the coupling is obtained also in this latter 
case. 

This paper is organized as follows. We describe the 
framework of our approach in sec. [2l We derive the main 
equation in sec. \3\ Its solution is presented in sec. [H We 
discuss its application in QCD and QED in sees. [SI and [SI 
respectively. The conclusions are stated at the end. 



2 General framework 

Let us begin with the gluon/photon Green's function: 

D{q,C) = Z{q) , (5) 

q-^ + te 

where ( is the gauge- fixing parameter and Z{q'^) is the 
renormalization coefficient. The running coupling, a(q'^), 
is proportional to Z{q'^). 

The central quantity which we adopt in the following 
discussion is the trace of the logarithm of D{q,(), which 
we denote by S^^. 

Hc(g)-TrlnI?(g,C). (6) 

Although this quantity will be shown to be essentially 
gauge invariant, we introduced the index C. to indicate 
that includes an unphysical polarization contribution. 

We define the logarithm in eqn. [6] by its Taylor series 
expansion, i.e.: 

In [A(/ + M)] ^ nnX + M - — + — -... , (7) 
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where A is scalar and M is a d x d matrix. / is the d- 
dimensional identity matrix. Applying this to eqn. [6l we 
obtain: 

This result is for the 4-dimensional vector boson, but is 
easily generalizable to arbitrary particles. Since the C de- 
pendence can be absorbed in the definition of Z{q^) or the 
scale of g^, is essentially gauge invariant. 

corresponds diagrammatically to the self-loop of 
the gauge boson, i.e., just a circle, without the phase-space 
integration. This assignment is natural because when we 
take its derivative, we obtain the trace of a propagator D 
attached to a 0-momentum vertex, —dD~^. Explicitly: 

d,S^{q)=Tv[-d^D-\q)D{q)\, (9) 

which follows directly from eqn. [SI We represent this rela- 
tionship diagrammatically as: 

(10) 

The thick solid line represents D{q, Q and the dashed line 
represents —dD~^[qX)- 

We would like to note at this point that the quantity 
inside the trace on the right-hand side of eqn. [9] has anal- 
ogous form to (q) introduced in ref . : 

A^{q) = d^G-\q) G{q). (11) 

Since we have: 

d^G{q) = -G{q)d,G-\q)G{q) = -G{q)A^{q), (12) 

the derivative of a propagator is equivalent to its multi- 
plication by —Af^{q) and is, by the Ward-Takahashi iden- 
tity, equivalent to the insertion of the 0-momentum pho- 
ton/gluon emission vertex. It is hence appropriate to rep- 
resent S'j(g) as a self-loop. 

In ref. T , (q) was introduced to reduce eqn. [3] to a 
first-order equation, namely: 

dMi) + (1 - 9)M<})Mi) = 0- (13) 

We reproduce this equation here for comparison with the 
photon/gluon Gribov equation, which we shall write in a 
similar form later. 

As for the physical interpretation of S'^, the logarithm 
of the propagator has the interpretation as the density of 
states |7]. Indeed, we see that the phase-difference: 

-lm[S^{q,)-S^{q,)], (14) 
TT 

by the Levinson theorem, is the number of states in that 
phase-space interval. However, as stated earlier, S'^ in- 
cludes an unphysical polarization contribution, and the 



number of states should in fact be 3/4 of this. This can be 
achieved by eliminating the scalar polarization component 
in eqn. [6l 

S(g)=TrP(g)lnC(g,C) = 31nf^^y (15) 

yq'^ +ie J 

Here, P is the transverse projection operator: 

P\a = gxa (16) 

q 

S{q) is gauge invariant. 

We may equally have chosen the derivative of S as the 
starting point of our discussion. Let us define: 

A^{qX)=d^D-\qX)D{qX). (17) 

We also define its transverse component: 

Al{q)=PA^{q,C)P, (18) 

which is gauge invariant. We can then define S as the 
indefinite integral of , i.e.: 

S{q) = -TT j' Al{q',Odq',. (19) 

Note that this is a line integral and not a phase-space 
volume integral which S{q) certainly is not. In general, 
because of the singularities in D{q^ C), the value of S de- 
pends on the path of integration in the complex plane. In 
particular, if we adopt a closed contour for the integration, 
we have an Aharonov-Bohm-type phase corresponding to 
the 'gauge field' TrA^(q). The value of the phase is 2i7r 
times the number of states enclosed by the contour. 

Using the above tools, the discussion of the Ward- 
Takahashi identity becomes relatively simple. To see this, 
let us write the 0-momentum limit of the bare 3-point 
vertex using the usual Feynman rules as: 

-rM,bare(<?, 9, 0) = 2gp.gA(T - 9A.9p<T - q^g^iX- (20) 

We have omitted the colour matrix f^^'^ for the sake 
of simplicity. This will be reintroduced later, in eqns. [35] 
and [Sol in the form of the colour factor Ga inherent in 
the beta-function coefficient bo. D~^, whose derivative we 
shall now compare with eqn. 1201 is given by: 

D-\q, C) - Z-\q') [-q'g^, + (1 - OqxQ.] • (21) 

Taking the bare propagator Dq^ , i.e., without considering 
the contribution due to the renormalization Z(q^), the 
derivative is given by: 

dD^\q, C) = ~2qf,gxa + (1 - C) (gAff^- + 9-3ma) • (22) 

The expressions [20| and [22| are related by the following 
tree-level identity: 

r^{q,q,0) = - (d^D-\qX) ~ C-^D-\q,C)) , (23) 
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exactly as found by Gribov 1 . The first term on the right- 
hand side gives the Ward-Takahashi identity, but the sec- 
ond is the compHcation arising from Slavnov- Taylor iden- 
tity. The latter term vanishes only for C = 0, but the 
transverse parts are the same. 

Now let us consider operations on E!{q). Since this 
is gauge invariant, we may choose the transverse gauge, 
=0. In this case, the insertion of a zero-momentum 
vertex anywhere in the all-order diagram for S'(g), i.e., 
including all vacuum-polarization contributions, is accom- 
panied by two transverse projectors, and therefore the 
non-transverse part of the vertex do not contribute. We 
hence seem to recover the Ward-Takahashi identity. 

This is not the end of the story, because even though 
the quantity D{q)r{q, q, 0)D{q) is transverse, the quantity 
lim^-i^o [Diq)dD^^{q)D{q)] is not necessarily so. The 
trick is to invoke the transversality of S'(g). For instance, 
we have: 

-Tr [D{q,C)dD~\qX)] - d^Sdq) (24) 
which differs from 



TT[D{q,C)r^iq,q,0)]^d^Siq), 



(25) 



but so long as we make use of E:{q), i.e., proceed by tak- 
ing the transverse component of the expression inside the 
trace in eqn. [24l the non-transverse component does not 
contribute and the two expressions become identical. In 
summary, the Ward-Takahashi identity: 



rn{q,q,0) 



\qXl 



(26) 



holds provided that we work with gauge-invariant and 
transverse quantities. For transversality, we look for the 
vanishing of the diagram when any gluon in it is assigned 
scalar polarization, ^^{k) — *■ k^. This holds for S{q) be- 
cause of its symmetry: the graph looks the same (and 
transverse) from anywhere in (the perturbative expansion 
of) the diagram. 

Unfortunately, in this discussion, we have lost gauge 
invariance: the transversality argument works only when 
we utilize a transverse gauge, C"""^ = with transverse 
S{q) and Al{q). 

A simpler working rule, which reproduces the same 
result provided that we are only interested in products of 
D and D^^ , is simply to work with the Feynman gauge 
and impose transversality at the end. To illustrate this 
point, A^((7, 1), which is written in the Feynman gauge, 
can be converted to Aj^{q) simply by multiplying on either 
side by the transverse projection operator P as can be 
easily verified: 



PA^{q,l) = A^iq,l)P 



(27) 



We have only discussed the tree-level relationship up 
to now, but the loop corrections decorate both sides of 
eqn. [21] in the same way, and so eqn. [21] is an all-order 
relation. 



3 Derivation of the equation 

In the Feynman gauge, the zeroth-order Gribov equation 
comes out immediately from eqn. [TJ 



d^D{q,l)KO (q^O). 



(28) 



In the rest of the discussion, we omit the second argument, 
C = 1, for simplicity. The small contribution to the right- 
hand side is due to the running of the coupling, so that it 
is 0{r), where P is the anomalous dimension defined by: 



d\nZ 
din (p- 



(29) 



This is related to the usual RGE beta function fi\s)j P — 
fi/a = d\na/dlnq^ . There is a difference of factor 1/2 
between our convention and that of ref . [1] . 

Eqn. [28] can be established in a more gauge invari- 
ant manner by using introduced in eqn. [1] Because 
of its gauge invariance, its derivative, — TrA^(q, C) is also 
gauge invariant. We are thus justified in using the Feyn- 
man gauge. Omitting the trace for simplicity, we can then 
write: 



dM<l) = d {dD-\q)D{q)) = -d {D-\q)dD{q)) 
Since: 

dD{q) = -D{q)dD-\q)D{q), 
and making use of eqn. 1281 we obtain: 



A^{q)A^{q)+0{P). 



(30) 
(31) 

(32) 



For the 0{P) term, since it is due to the running of 
the coupling, we need to analyze the polarization operator 
n{q^) = — i7pp((7)/3(7^. Let us see how this enters our 
framework. 

The expansion of S'(g) yields the following: 



(33) 

The thin lines represent the bare propagator Do{q), and 
the blobs represent the insertion of the polarization op- 
erator. The factors 1/n are the symmetry factors for the 
circular symmetry. The contribution of the polarization 
operator, i.e, the sum of the contributions from all but 
the first term of eqn. I33i is proportional to: 



7T(g2) 



= -ln(l-77(g2)). (34) 



Including the first term of eqn. 

(1- 



we thus obtain: 



S((j) = 31n 



i7(g2))- 



■ te 



(35) 
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Comparing with eqn. 1151 we obtain: 



(36) 



which is as expected and consistent with the usual renor- 
malization considerations. 

In order to find the sub-leading term in eqn. I32( we 
take the derivative of eqn. 1331 Taking the first derivative, 
we have a series of diagrams, approximately half of which 
involving the derivatives of the line and the rest involv- 
ing the derivatives of the blob. The derivatives of the line 
are already taken account of in the leading-order term of 
eqn. [32l and so we are left with the following series of 
diagrams: 



0-0*0 



(37) 



The combinatorial factor at each order cancels with the 
symmetry factor present in eqn. 1331 

When we take the next derivative, we again have some 
derivatives of the line, and these are already included in 
eqn. 1321 In addition, we have cases in which the two deriva- 
tives act on separate blobs, such as: 



this operation by 9^, where the subscript R stands for 
renormalization. With this understanding, we may write 
the equation in terms of n{q^), as: 



1 



Tr [dMq) - M<i)Mq)] = o^dln{q^ 



(40) 



The left-hand side of this equation is not gauge invariant, 
and the Feynman gauge is implied. It is a trivial matter to 
convert this into a transverse expression, just by substi- 
tuting A^{q) for Ap(g) and averaging factor 1/3 for 1/4. 

It is interesting to explicitly evaluate ad^n{q^) and 
demonstrate that this does not equal to the left-hand side 
of eqn. |40l To do so, let us first write out the explicit form 
of A^{q). For general choice of C, we have: 

A^{qX) ^ d,D-\qX)D{qX) 
= %(l-r).gA.+ 

^iC-i)P,. + ^ii-C')P,x. (41) 

We see that the gauge dependent terms, which are propor- 
tional to q\ and qa-, vanish when either we take the trace 
of A^, or in the Feynman gauge. In the Feynman gauge, 
we have: 



A^(g) = %(l-r)gAa 



(42) 




(38) 



These give rise to the renormalization of eqn. 1321 and 
therefore are already included. With this consideration, 
the remaining second-derivative terms are: 

0-0*0* o 

(39) 

As before, the thick solid line represents the renormalized 
propagator. As for the blob, it is essentially Uxcr- However, 
evaluating d^TIxa will not yield a sensible answer, because 
there are cancellations between the propagators ~ 1/q^ 
and JTact ~ 9^ which make the theory renormalizable and 
make eqn. [32] valid. 

This ~ q^ factor, in the Feynman gauge and at the, 
dressed if necessary, one-loop level (which is all that is 
needed), arises from the linear momentum dependence of 
the external vertices, i.e., the vertices which are connected 
with the external propagators. Thus we need an operation 
which extracts only the contributions to d^IIx^ that arises 
from the derivatives of the (dressed) internal lines and 
not the oc q linear external-momentum dependence of the 
vertices attached to external propagators. Let us denote 



It is a simple matter to work out A^^{q) and d^A^{q). 



These are given by: 



A,{q)A,{q) = {l-rf^., 



and 



r-(i-r) 



(43) 



(44) 



respectively. Here, F refers to the derivative by In (7^. We 
may also show that exactly the same relations hold for 
A^{q), provided that we replace g\a on the right-hand 
side by Pa<t- 

From eqns. [43land[44l we obtain: 



hYr[dMq) - A,{q)A,{q)]^ ^ 



-r + r- 



On the other hand, ad^n{q^) yields, by using eqn. 



(45) 



ad'n{q') = -Z{q')d'Z-\q') = - 



r + r -r' 



(46) 

i.e., similar to eqn. 35] but with a different sign for F. 
Therefore d\n{q^) and d^n{q^) differ. 

Now turning our attention to 9^77(g^), and temporar- 
ily adopting the language of perturbation theory with bare 
propagators, n{q^) has an expansion containing, in gen- 
eral, many gluon lines. 

Following Gribov's argument mentioned in the intro- 
duction, the logarithmically enhanced contribution to the 
polarization operator comes from regions in the phase 
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space which involve large hierarchy of internal momenta. 
The greatest contribution to d\n{q^) comes from terms 
in which both derivatives act on the same line, and so we 
may apply eqn. [Tl which removes a momentum integra- 
tion. 

This operation converts the amplitude into that of the 
emission of two zero-momentum gluons, so that the lead- 
ing contribution will be given by: 

dln{q^) cx Tr [A^{q)A^{q)] + 0{r), (47) 

including all the renormalization correction to the vertices 
and propagators. 

At the one-loop order, we have, independently of the 
gauge parameter: 



1 



-47r2i -I- ie){{q ~ k)^ + ie) 



, (48) 



where the subscript in 77i refers to the perturbative order. 
bo is the first beta function coefficient, which we define to 
be positive for asymptotically-free theories. Applying 
to this expression, we obtain: 



dln.iq) = d^n,{q) - -46o^— -. (49) 

q^ + le 

Comparing with the form of eqn. l43[ this fixes the constant 
of proportionality in eqn. [47lto be —bo /A. 

Now let us work out the same constant of proportional- 
ity in an all-order approach based on dressed propagators. 
Inserting (some of the) Dyson-Schwinger type all-order 
corrections, as shown in fig. [1] into eqn. 1481 we have: 

b f k 

nsuiq) = -x^Tr / ——rx.h.reD{k)nD{q^k). (50) 
As in eqn. I20| /bare represents the bare vertex. 



q — k 



Fig. 1. The Dyson-Schwinger- type correction to the one- loop 
gluonic vacuum-polarization operator. The thick lines and the 
vertex to the right are renormalized. The thin line and the 
vertex to the left are bare. 



As it stands, eqn. [50] is problematic in the sense that 
it is not explicitly gauge invariant, and if it is, we should 
also include the contribution due to the ghost. 

Since we would like to make use of the Ward-Takahashi 
identity of eqn. 1261 it is necessary to work with transverse 
quantities. As discussed in sec. [H a simple way to do this 
is to work in the Feynman gauge and impose transver- 
sality in the end. It is in fact not even necessary to rein- 
troduce transversality at the end, but as a penalty, the 
spin-averaging factor will be replaced by 1/4. 



We may thus replace Fx by —d\D^^ (Feynman gauge), 
with the understanding that the phase-space region giving 
rise to large logarithms has q — k ^ q, k and hence the 
zero-momentum emission vertex is a good approximation 
to the full vertex. In this case, the simplest choice of scale 
would be D~^{k). The integrand in eqn. [50l then becomes: 



dxD^'{k)D{k)dxD-Hk)D{q - k). 



(51) 



When we apply dj^ to eqn. [501 the leading renormal- 
ization contribution comes from the double derivative, d^, 
of D{q — k). The application of to D{q — k) yields a 
delta function in the approximation that the coupling is 
constant. If not, we have a correction term proportional 
to r. The first derivative of D{q — k) yields: 



2Z{1 - r)q^ 
{{q - k)'^ + ie)' 



(52) 



Hence d^D(q — k) is approximately multiplied by factor 
—Z{\ — r) as compared with eqn. [TJ Factors of Z cancel 
in the expression for djiUgD, which now reads: 



Mi—Qxr [dxD,\q)Do{q)dxD-\q)D{q)] . 



(53) 



Let us make the scale choice as r{q'^). The bare propagator 
Da{q) and the renormalized propagator D{q) differ only 
by the factor Z. Similarly, the difference between dxD^^ 
and dxD-'^ is, by eqn. [42l Z'^il - F). Thus we have: 

dln^^iq) ^ -^Tr [dxD-\q)D{q)dxD-\q)D{q)\ . 

We have now shown that by applying Dyson-Schwinger- 
type corrections to the vacuum-polarization graph, we are 
able to reproduce the form which is expected by the log- 
arithmic enhancement argument of Gribov which implies 
eqn. [ITj without the need of such algebraic manipulations. 
Although this correspondence between the two approaches 
may seem intuitive and natural, we are not sure about how 
one would go about formulating such a correspondence in 
the case of the fermionic Gribov equation, eqn. [3] 

In any case, together with eqn.l40|, we have now estab- 
lished the following equation: 



dMD = (1 - &oa) A^{q)A^{q) + o{r^). 



(55) 



According to eqn. 1321 the term in eqn. [55] that is pro- 
portional to 60^: which corrects eqn. [32l is 0{r). This is 
true, since —bo is the first expansion coefficient of the beta 
function. We have: 



boa^r + oir^). 



(56) 



Taking the first term, we obtain a rather compact expres- 
sion with no parameter dependence: 

dMi) = (1 + nq')) Mi)M<i) + o{r^), (57) 

from which we expect bo to reappear as a constant of in- 
tegration. We shall see in the paragraph following eqn. [76] 
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that this reproduces the magnitude of the next-order co- 
efficient, b' , of the beta function expansion. Therefore the 
error in eqn. [57] is practically 0{r^). It is actually dif- 
ficult to modify eqn. [56] without introducing unphysical 
fixed points. This point will be discussed further in the 
paragraph leading up to eqn. [60] 



4 Solution of the equation 

Let us now solve our Gribov equation, eqn. 1571 

Substituting eqns. [43l and [44] into eqn. [57] we obtain: 

r-(i-r) + (i + r)(i-r)2 ^ r-r^{i-r) = o. (58) 



That is, 



d{l3/a) 



(59) 



r = P/a is positive for QED and negative for QCD. 
The equation is nominally not applicable to the case of 
QED beyond the one-loop order, or to the vacuum polar- 
ization due to the quark loop in QCD. On the other hand, 
we believe that the formalism is still useful in describing 
the UV behaviour of QED because of the fixed points in- 
herent in eqn. [59] Note that the zeros of the right-hand 
side of eqn. [59] represent fixed points. The F = fixed 
point corresponds to the trivial vacuum whereas F ~ I 
corresponds to the QED UV fixed point. In this limit, the 
running of the QED coupling cancels the propagator 
factor, and so the photon decouples. 

Even though it is not obvious that the photon should 
necessarily decouple in the UV limit, a limit in which the 
photon decouples is almost necessarily a fixed point of 
the theory, since in this limit there is no longer photon 
propagation and hence no further evolution of the photon 
propagator. In this regard, the F = 1 fixed point of eqn. 1591 
seems physical. 

On the other hand, the equation will certainly break 
down when discussing effects due to fermion masses or 
fermionic super-critical state formation. Related to this 
point, we do not expect eqn. [55] to be valid when the ex- 
pansion parameter, i.e., in this case, F, is large. 

In the case of UV QED, we think that the evolution 
is barely permissible because of the presence of the phys- 
ical fixed point, but in the case of IR QCD, as there is no 
fixed point for negative 7^, 7^ diverges at, as we shall show, 
^QCD ■ In this case, we do not believe that the equation is 
quantitatively correct. However, the behaviour of the solu- 
tion is, we believe, nevertheless physical and, in any case, 
Gribov's super-critical state formation occurs before this 
singularity. A measure of the quantitative accuracy of the 
equation is provided by the comparison with the pertur- 
bative beta-function expansion, and this will be presented 
in the next section. 

We would like to note, to avoid confusion, that we also 
expect the F — 1 fixed point of eqn.[59]to be correct in the 
case of QCD because of the decoupling behaviour which 
it represents. However, the limit indicated by this fixed 
point does not arise in asymptotically-free theories. 



Before proceeding, we note that, had we chosen a dif- 
ferent expression for eqn. 1561 we would, in general, end 
up with extra fixed points with unphysical power-like be- 
haviour of the coupling. Because of this, it is difficult to 
introduce a simple higher-order modification to eqn. [55] 
without affecting its physical behaviour. On the other 
hand, if there arises a need to create a toy model for the 
running coupling with some particular power-like fixed- 
point behaviour, it is easy to artificially modify eqn. 1561 to 
serve this purpose. For instance, one may like to introduce 
an which is finite in the space-like region. A possibility 
would then be the substitution: 



boa — F 



F^ 
1~F' 



(60) 



In this case, corresponding to eqn. 1591 we obtain: 

F ^ F^{1- F^), (61) 



which is self-dual under Inq^ 



In q'^ and F 



-F. In 



the low-energy limit of QCD, this gives us a a oc 
behaviour, i.e., a simple toy model for the long-distance 
linear potential. 

Resuming our discussion of eqn. 1591 let us first confirm 
that it leads to the ordinary result for the running cou- 
pling when 7^ is small. We would like to calculate the evo- 
lution in the space-like region, i.e., for positive = —q^. 
Omitting the sub-leading term, we obtain: 



P/a = [- InQ^ + const. 



(62) 



The constant on the right-hand side is In . The left-hand 
side is the logarithmic derivative of In a, so we obtain: 



lna = J [\n{A^/Q^)] ^ d\nQ^ 
= - ln(ln(ylVQ^)) + const. 



(63) 



This second integration constant gives the leading-order 
beta- function coefficient bo, so we finally obtain: 



am 



1 



boHQ^A^)- 



(64) 



As before, bo is, in our convention, positive for QCD and 
negative for QED. 

More generally, the integration of eqn. [59] yields: 



ln(ylVQ' 



F- 



In 



1 



(65) 



The modulus is a shorthand for writing 7^^^ — 1 for QED 
and 1 — 7^~^ for QCD. Even though it is not possible to 
invert this equation, there is a trick to integrate it. We 
have: 



log a 



FdlnQ^ 



F'-^dF^ [LdF 



dF 



F 



(66) 



Then by the use of eqn. [59] we obtain a relation between 
a and 7^: 



log a = 



dF 



F{l~F) 



(67) 
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6q ^ is the constant of integration. Thus we have: 

(boa)-' = l-r-\ (68) 

which implies the following beta-function expansion: 

I3 = ar = -boa^ (1 + boa + bla^ + •••). (69) 

We shall compare this with the QCD beta-function coef- 
ficients in the paragraph following eqn. 1761 We now sub- 
stitute eqn. [B5]in eqn. [^5] to obtain: 

In(QV^e) = (M)"'+ln|^>o|a, (70) 

or: 

QVyl2 = |6o|aexp(l/6oa). (71) 

We have defined = Ae~^^^. In the following, let us 
discuss the properties of this solution. 



The QCD coupling 



1 




Real part 

Imaginary part ■ 

1/log(dV) 
1/[3+log(Q^/A^)] 





10"" 10"^ 10° 10^ lO" 



Fig. 2. The QCD running coupling calculated using the Gri- 
bov equation, compared against the leading-order perturbative 
behaviour, and the same shifted by yl^ ^ yl^e"^ 



5 The case of QCD 

Eqn. [55] and its solution, eqns. [7D] and [7T1 are general 
both to QED and QCD. However, eqn. [59] has no zeros 
for negative F, and so its solution becomes non- analytical 
at = yl^ = ^QCD- this point, boa reaches 1. 

To see the behaviour of eqn. [71] below Aqqd , let us 
write: 

boa = ae"''. (72) 
Then by taking the imaginary part of eqn. 1711 we obtain 

a<f> = sin (j). (73) 

At = yl^, a = 1 and = 0. Below Q'^ = A^ , depending 
on whether we move the singularity above or below the 
real axis, becomes positive or negative. Adopting posi- 
tive (/), a gradually decreases, and by the form of eqn. [TS] 
4> gradually increases. In the limit — s- 0, a vanishes, 
and we end up with (j) = -k, i.e., a becomes negative. This 
rotation of phase implies that there is one gluonic super- 
critical state between = and = A, with negative 
mass 13 . 

Taking the real part of eqn. [TT] we obtain the other 
constraint: 



Al 



■ exp 



tan ( 



(74) 



We can now make a plot of boa against Q"^ / A for both real 
and imaginary regions. In the real region, we use eqn.lTllto 
evaluate Q'^ / A as a function of boa, whereas in the imag- 
inary region, we make use of eqns. [73l and (74] to evaluate 
a and A^^/Q"^ as a function of 0. Fig. [2] shows the plot 
obtained in this way. 

Above — A'^, we see that there is a modification 
to the one-loop perturbative evolution, which persists up 
to considerably high energy, but most of the modification 
can be absorbed by a shift in ^qcd , found numerically to 
be about ~ 7. This implies that the measured Aqcb 



is about three times smaller than the true Aqcd that is 
obtained in the high-energy limit. 

To study this effect, let us consider the iterative inver- 
sion of eqn. [TO] 



(boa) ^ = log 



Q 



log(...) 



(75) 

At the first level of truncation we obtain A^ — > A^boao 
where ao is the value of a at some relevant scale, and so 
on. 

In fact, such a shift in /Iqcd is known to be present 
already at the perturbative level. One expression for as 
at the two-loop order reads [5]: 



b' In 



1 + b'asiQ^) 



boHQ^/A"), (76) 



where b' is the ratio of the first and second beta-function 
coefficients: 

/9porturbativc = -boa^{l + b' a + b"a^ + •••)■ (77) 

We see that eqn. [76] has almost the same form as eqn. [70] 
Since the presence of an extra a inside the denominator 
of the second term is a higher-order effect, and so is the 
choice of A, the two equations differ only by the difference 
between bo and b', as can be inferred from eqn. 1691 In the 
real-world QCD, this difference is given by: 

b' _ 6(153- 19?!/) 
b^ ~ (33-2n/)2 

for Uf = 3. On the other hand, the large- A^c or n/ = 0, 
limit of the same quantity is 0.843 • • •. In either case, it is 
reasonably close to unity. 

To analyze the higher-order contributions [^lOj. let 
us introduce an alternative and common notation for the 
beta-function coefficients: 



0.790- ••, 



(78) 



P/a 



/3„_i(a/4^)". 



(79) 
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Uf ^3 




Nc = 00 


P2 = {4nfbob" 
133 = {4n)V 


9 

64 
3863/6 
12090.4 


11 

102 
2859/2 
29243.0 


102(iVc/3)^ 
2859(iVc/3)V2 
25554.8(iVc/3)^ 


Pi/Po 
132/ f3i 
P3/P2 


7.11 
10.06 
18.78 


9.27 
14.00 
20.47 


9.27(Afc/3) 
14.00(iVc/3) 
17.89(iVc/3) 


{P2/Pl)/{Pl/Po) 
{P3/P2)/W2/Pl) 


1.415 
1.867 


1.510 
1.462 


1.510 
1.277 



Table 1. The beta function coefficients and their ratios. The 
three-loop [5] and four-loop [10| contributions are calculated 
in the MS scheme. 



These coefficients are shown in tab. [l]for the three cases: 
three light-quark flavours, zero flavours, and for Nc — 00 
which implies zero flavours. The zero-flavour case starts to 
differ from the Nc = 00 case only at the four-loop order, 
where the non-planar contributions which are absent up 
to the three-loop order arise [10] . (32 and are dependent 
on the renormalization scheme, and the numbers quoted 
in tab. [1] correspond to the MS scheme. The choice of the 
renormalization scheme will affect our discussion here, at 
least in principle. 

One notices that even though the ratio between two 
successive coefficients, /3„//3„_i, deviates away from /?o at 
higher orders, it seems to tend to a constant, when the 
number of flavours is zero. Of course, it is dangerous to 
draw any conclusions through knowing only these four co- 
efficients, and our knowledge about the higher order co- 
efficients is limited, but let us proceed with this tentative 
discussion for now. 

To quantify this statement about the ratio of two suc- 
cessive coefficients tending to a constant, we have also 
tabulated the ratio: 



(/3„+i//3„)/(/3„//3„-i) = /3„+i/3„_i//32, 



(80) 



in tab. [TJ For the beta-function expansion to have a finite 
radius of convergence, by d'Alembert's ratio test, it is nec- 
essary, though not sufficient, for this ratio to tend to one. 
The point at which the beta-function expansion first di- 
verges is the point at which super-critical behaviour arises. 
We expect the ratio to be positive, since its being negative 
would imply super-criticality for negative as, that is, for 
repulsive strong interaction, and this is unphysical. 

The values of (3n+i/(3n being different from (3q implies 
that eqn. [69] is obviously not literally correct, but is an 
approximation in the sense that up to the three-loop or- 
der, these ratios are actually close to /3o, even with finite 
(small) number of ffavours. This gives partial assurance 
about the validity of our approach. 

At higher order than three loops, in our opinion, the 
validity of the approach lies not in the numerical accuracy 
of the Gribov-equation evolution but rather in the physics 
it describes, namely super-criticality at finite as- 

As for the formal presence of the 7^ = 1 fixed point 
in eqn. 1591 this means that in the limit of large as, when 



the /3-function is analytically continued beyond its formal 
radius of convergence, we expect P/a to tend to 1. 

Let us return to the discussion of the branch-point 
singularity mentioned at the beginning of this section. The 
formation of the gluonic super-critical state, indicated by 
this singularity, would make the vacuum unstable, since 
there is energy gain in, for example, the pair production 
of such states from the vacuum. As far as we can see, there 
is no way to stop or saturate this decay, either statistical 
or dynamical, and so pure QCD has to be, in our opinion, 
inconsistent. 

However, the situation becomes different when we have 
light quarks. According to Gribov, the bound states of a 
pair of light quarks becomes super-critical when the cou- 
pling exceeds the critical coupling ac given by: 



(81) 



This is smaller than the gluonic branch point which occurs 
at boa = 1: 



-Cf 



16 



fcna— 1 



33 - 2nf 



0.592- 



(82) 



As before, we have taken Uf =3. We note that, poten- 
tially, the effective coupling which appears in the Gribov 
equation of eqn. [3] should be corrected by (1 — r'), which 
arises from eqn. [55] This is a small effect for the value 
of as given by eqn. [HU but it reduces the critical cou- 
pling. In addition, the actual values of the beta function 
coefficients, or the ratios of them, tabulated in tab. [Jim- 
ply that the gluonic super-criticality occurs at much lower 
values of as, probably a half or so, than is implied by 
60a = 1. Explicitly, by the ratio test for the convergence 
of the beta function expansion, giuon super-criticality oc- 
curs at as — af'" given by: 



/An — lim 



/3. 



N 



N- 



P. 



N+1 



(83) 



It is not clear per se whether this limit is actually finite, or 
that it is renormalization-scheme independent. However, 
our work suggests that it is at most finite and is proba- 
bly non-zero. If it is zero, the implication would be that 
there is some form of super-criticality at any value of the 
coupling, and this appears unlikely to us. 

Returning to the discussion of the presence of the light 
quarks, since it is a non-trivial matter to consider the 
coupled evolution of the gluon and quark Green's func- 
tions, let us assume that Gribov's argument leading to 
the formation of light-quark super-critical states is essen- 
tially unmodified by the running of ag. Obviously this is 
questionable if the singularity at — remains, but 
we believe that this is not the case. 

Firstly, due the formation of the light-quark super- 
critical state, the quark Green's function becomes com- 
plex, and this corresponds to the decay of the light quarks 
by emitting the super-critical state, or more strictly the 
decay of the vacuum. When this occurs, the gluon Green's 
function also becomes complex, indicating the decay of 
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the gluon by emitting the super-critical state through the 
quarks. This would move the singularity off the real axis, 
and the evolution of ag will continue to = without 
any further singularities. 



( 




•2 

» • - 




1 



Fig. 3. The singularities of the gluon Green's function. The 
branch-point at Q'^ = is labelled 0, that at A? is labelled 
1. 2 corresponds to 1 moved off axis due to the decay into 
light-quark super-critical states via light quarks. 

This situation is shown in fig. [S] In one-loop perturba- 
tive QCD, there is a branch-point at + ie = 0, as well 
as a simple pole at = —q^ = ^qcj). However, the non- 
perturbative effects inherent in eqn.[57]tame this and make 
it a branch point, corresponding to gluonic super-critical 
singularity. The decay into light-quark super-critical states 
move this singularity off the real axis. 

This is when the negative-energy super-critical states 
are not completely filled. When they are filled, these de- 
cays become forbidden, and so as becomes real again. 
When this occurs, the singularity near — vlq^p should 
disappear, but there will be a new branch cut starting at 
q^ — ie — which, as in ref. [T] , arises as a result of what 
is best described as the reorganization of the Dirac sea. 

To describe the evolution of the Green's functions with 
the new Dirac sea, according to the approach adopted in 
ref. [T], it is sufficient to modify the evolution equations 
by including the contribution of the Goldstone boson (the 
pion) . Since the Goldstone boson only couples to the gluon 
through quarks, the inclusion of the Goldstone boson con- 
tribution would be through the light-quark loops in the 
vacuum polarization operator. 

As stated above, the coupled evolution of quark and 
gluon Green's function is a non-trivial matter, but so long 
as the modification is local in the momentum space as in- 
dicated by ref. [T], the IR behaviour of asiQ"^) is governed 
by eqn. [59] Even if the locality does not hold, in general, 
we expect the quark contribution to the running of the 
QCD coupling to be small. Here it causes a large effect 
only indirectly, through (the forbidding of) the decay of 
the gluon. 

From the locality of the correction due to the Gold- 
stone boson and the analyticity of ces{Q^), we may deduce 
the following. 

If 7^ remains negative below — ^qcq, a branch- 
point singularity in as{Q^) again develops at finite Q^. 
Since this would contradict analyticity, F must change 
the sign. It must also remain below F — I, because other- 



wise it again develops a singularity when F diverges. As 
mentioned at the beginning of sec. [H we expect that the 
F — 1 fixed point is correct even in asymptotically-free 
theories in which this UV fixed-point limit does not arise. 
We then have an IR-free theory in the IR limit. Fig. 2] 
illustrates this behaviour. 

as(Q') 



chiral- 

symmetry 

breaking 




Fig. 4. The behaviour of the QCD coupling at low energy. 

Since, as we mentioned above, we do not expect the 
quark (dynamical) contribution to the running of the QCD 
coupling to be large, the evolution in the IR-free region 
is almost entirely due to the gluons. Thus it may seem 
strange that the gluons screen the coupling here rather 
than yield the usual anti-screening behaviour. However, 
this is as expected, because in the low-energy effective 
theory, as a result of the forbidden decay to the super- 
critical states (or in other words, the decay into the new 
states created by Dirac-sea reorganization), there is a new 
cut with the branch point sX q^ — ie = 0. This means that 
the Wick rotation has to be performed in the opposite 
direction to the usual one, and so the direction of the 
running is inverted. 

This scenario is as discussed by Gribov [ll[3] . The new 
cuts, which apparently violate causality, correspond to 
the presence of the positive-energy super-critical states of 
negative-energy quarks. Causality is not violated, but it 
appears as if positive-energy quarks are travelling back- 
wards in time. 

An IR-free gluon is decoupled at = 0. Thus it 
cannot be found as a free particle, though bound states 
with finite radii can contain it. The masses of these bound 
states would be of the order Aqcb purely by dimensional 
considerations. Hence the energy scale, or the energy gap, 
for gluon decoupling is ©(/Iqcd)- This is soft confinement, 
meaning that the gluons are bound together only by finite- 
distance dynamics. 



6 The case of QED 

Let us now turn our attention to the case of the running 
couphng in QED. 

We show the plot of eqn. [7T] for negative bo in fig. O 
This plot is generated by evaluating Q'^/ as a function 
of boa 
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The QED coupling 



10° 



10=' 



10" 



10^ 



10' 



10' 



10" 



10-' 



10"' 



lo--^ 



10"' 




Q'/A" 



Fig. 5. The QED running couphng calculated using the Gri- 
bov equation, compared against the leading-order perturbative 
behaviour and the oc behaviour. 



The high-energy behaviour is governed by the F — \ 
fixed point of eqn. [59l The large-couphng limit of eqn. [71] 
is given by: 

lim a{Q^) = 



\ho\Al' 



(84) 



This corresponds to the straight line shown in fig. [5l We 
also show the ordinary one-loop perturbative result in the 
same plot. 

As in QCD, we expect that the value of A is highly 
sensitive to various factors, including the higher-order cor- 
rections. Therefore we do not trust the constant of pro- 
portionality in eqn. [84] to be accurate. On the other hand, 
we believe that the oc behaviour to be correct since, 
as explained earlier, the photon decouples from the theory 
only in this case. 

This oc behaviour for the UV couphng yields a 
current-current contact interaction, similar to the NJLVL 
model [6]. The couphng constant, M^jlvli by eqn. [HH is 
given by ATrboA^. The interaction Lagrangian density is, 
for space-like exchange: 



1 



NJLVL 



(85) 



Ref. [6] reports the formation of a massless Goldstone 
mode. This is consistent with Gribov's EWSB mechanism 
[S] based on top quark condensation [TTl due to the U{1)y 
Landau pole. There are other bosonic modes such as the 
Higgs boson and the massive axial vector boson reported 
in ref. [6j. Although the photon decouples at = yl^, 
these bosons remain physical even above this scal^. 



^ This behaviour of QED at high energy is similar to that in 
media with large screening. When either the coupling a or the 
response of the medium |&o| is large, the interaction becomes 
point-like, i.e., the contact-term interaction. The oscillation of 
the medium gives rise to states like the charge density wave 
in solids and the Higgs and Goldstone bosons in high-energy 
QED. 



Our results are supported by a numerical study of 
high-energy QED JJj, which yielded exactly the same 
conclusions, namely that the photon decouples, yielding 
a contact interaction which gives rise to chiral-symmetry 
breaking. 

We believe that the effect of gravity does not spoil 
the applicability of this EWSB mechanism, because the 
relevant gravitational coupling remains small above the 
Planck scale |13j . 

There are no fixed points for time-like running, and 
therefore eqn. 1641 remains a valid description for negative 
— —q^. In this case, the real part of the couphng 
changes its sign at = A? , indicating that the photon 
becomes a ghost, or unphysical, above this value of . 

Let us now compare our results with eqn. [4| due to 
Gribov. As mentioned in the introduction, the strong- 
coupling limit of this equation is given by: 



(86) 



where ^ = InQ^/yl^. If we omit the g ^ term, the general 
solution of this equation is: 



9 



(87) 



Ci and 6*2 are the constants of integration. The coupling 
asymptotically tends to C^^ at high energy. 

With the inclusion of the term, this is no longer 
constant, and we obtain — > C^^j or: 

g ^ IniQ^A'), (88) 

which is the behaviour obtained by Gribov [4]. 

Let us consider the modification to eqn. [551 due to the 
running of the photon propagator, viz. eqn. 1521 



9cs = .9(1 - r)- 



(89) 



We make this substitution because the g on the right-hand 
side of eqns. [D and [5S| arises from the solution of eqn. [31 
By exactly the same argument as that in the paragraph 
following eqn. 1521 there is an extra renormalization effect 
which multiplies g hy {1 — F) . 

Since the left-hand side of eqn. [86] is derived by manip- 
ulating the derivatives of n{q^), it is unaffected by this 
effect. We then have: 



d\l 2 

5? + ^5ejr-?(T^^ 



We have assumed that the strong-coupling approximation 
is valid in the sense that gcs is greater than approximately 
1. An asymptotic solution, this time, is: 



F^l- V^A^, 



9 



(91) 



This choice of F makes g{l — F)^ — > 2 asymptotically 
constant, and so g / A^ — solves eqn. [TO] 

Hence we believe that our results are consistent with 
that of Gribov, provided that one takes into account the 
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effect due to photon renormalization inside the vacuum- 
polarization operator. 

Before concluding this section, we would like to men- 
tion one property of eqn. 1711 which does not seem partic- 
ularly useful, but we think is worth mentioning. 

Obviously eqn. [71] cannot be inverted to obtain a as 
some elementary function of Q^, but certain moments of 
it can be evaluated in a closed form and are finite. Let 
X = Q'^/Al and a = \bo\ a. Then: 

/•OO POO 

/ a;"a-'"da= / e-"/"a"-'"da. (92) 
Jo Jo 

Then by the definition of the Eulcr Gamma function, this 
becomes: 

ni+"-'"r(m-n-l). (93) 

The moments of x" under a may not be useful at all, but 
they are related to the moments of a^™ under x, which 
may seem slightly more useful. 



7 Conclusions 

We derived a local Gribov equation for the gluon/photon 
Green's function D{q), and solved it for both QCD and 
QED. 

Our derivation is based on taking the second deriva- 
tive of Ti'ln D{q). We separated out the parts due to the 
running coupling from the parts due to the renormal- 
ization of the leading, two-gluon(photon)-insertion term. 
Both using Gribov's logarithmic-enhancement argument 
and using a Dyson-Schwinger-type expression, the part 
due to the running coupling is shown also to be of the 
form which corresponds to the double emission of zero- 
momentum gluons/photons, and so the Gribov equation 
can be written down in a compact form. 

The Gribov equation gives an equation for the running 
of the coupling in a closed form. We obtained the solution 
of this equation in an analytical form, for both QCD and 
QED. Although we do not expect our equation to be fully 
applicable to QED, we argued that it has sensible UV 
behaviour. 

In the case of QCD, we obtain an as which is finite, 
but has a branch-point singularity at = Aqq^^. We 
interpret this as being due to the formation of gluonic 
super-critical states, which makes the vacuum unstable. 
However, adopting Gribov's scenario, the singularity is 
moved off the axis due to the decay of the gluon into 
light-quark super-critical states. Furthermore, with the re- 
organization of the (Dirac-sea) vacuum, these decays be- 
come forbidden, giving rise to a coupling which only has 
singularities along the time- like axis. However, there are 
singularities on both +i£ and —ie sides of the time-like 
axis, giving rise to, in the low-energy limit, a coupling 
which is IR-free. The gluon is then softly confined. 

In the case of QED, the Gribov equation has both 
an IR fixed point given by -T = and UV fixed point 
given hy r = 1. The latter gives a UV oc limiting 
behaviour for the coupling. The photon decouples from 



the theory. The high energy limit of QED is then given by 
the contact interaction. This supports Gribov's scenario 
of EWSB by top-quark condensation due to the strong 
U{1)y interaction near the would-be Landau pole. 

We believe that our formalism can be applied also to 
the problem of scale generation in gravity [13] . 
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